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A tunable and balanced heat interferometer is proposed and analyzed. The device consists of two supercon- 
ductors linked together to form a double- loop interrupted by three Josephson junctions coupled in parallel. 
Both superconductors are held at different temperatures allowing the heat currents flowing through the struc- 
ture to interfere. As we show here, thermal transport is coherently modulated through the application of 
a magnetic flux. Furthermore, such modulation can be tailored at will through the application of an extra 
control flux. In addition we show that, provided a proper choice of the system parameters, a fully balanced 
interferometer is obtained. The latter means that the phase-coherent part of heat current can be controlled 
to the extent of being fully suppressed. Such a device allows for a versatile operation appearing, therefore, 
as an attractive key to the onset of low-temperature coherent caloritronic circuits. 



Manipulation of heat currents in mesoscopic circuits 
represents a youthful field of research.lH^The role of dissi- 
pative phenomena at the nanoscale and the benefits that 
can be derived from them are becoming to be under- 
stood. Starting from cooling applications and fine tem- 
perature tuning in cryogenic detectors,^ or superconduct- 
ing circuits for quantum computing,^ and ending with 
the emergence of caloritronic circuit s.^K^^ Towards this 
end, a magnetic flux- controllable superconducting heat 
interferometer was recently theoretic ally conceived^^ and 
subsequently realized experimentally.^ This achievement 
served to show, on the one hand, how quantum coher- 
ence between two weak-linked superconducting conden- 
sates extends also to dissipative observables such as heat 
current. On the other, this device might constitute the 
building block for the implementation of superconduct- 
ing hybrid coherent caloritronic circuits like, for instance, 
thermal modulators, heat transistors and splitters, etc. 
In this Letter we propose a step forward this goal by envi- 
sioning and theoretically analyzing a fully-balanced heat 
interferometer. This device provides an enhanced control 
over the flux-to-heat current transfer function therefore 
enabling the choice of convenient operation points for 
different applications. We show as well how, under ap- 
propriate conditions, the phase-dependent component of 
heat current can be maximized or completely annihilated. 

Our thermal circuit consists of two superconductors Si 
and S2, weak linked forming a double- loop interrupted 
by three parallel Josephson junctions (see Fig. [T]). Let 
us denote Ri the normal-state resistance of junction z, 
and (^i the macroscopic phase difference across junction 
i. This structure behaves as a conventional supercon- 
ducting quantum interference device (SQUID) pierced 
by a magnetic flux <l>i in which one of the junctions 
has been replaced by a DC SQUID. The characteris- 
tics of this second "junction" can be tuned thanks to 
the application of a control magnetic flux ^2- The sys- 
tem is temperature biased by setting the temperature in 




FIG. 1. Our superconducting circuit consists of three par- 
allel Josephson junctions (characterized by their Josephson 
currents ia, ih and ic, respectively) that define a double- loop 
heat interferometer. The temperature in Si is risen up to 
> which yields a steady-state heat current Jtot flow- 
ing from Si to S2. Moreover, the voltage drop across the 
device is set to zero. The main loop is pierced by a mag- 
netic flux $1 whereas a control flux $2 is applied through the 
second loop enabling enhanced performance. 



Si to be Tsi > Ts2^ Ts2 being the temperature in S2. 
Furthermore, the voltage drop across the whole struc- 
ture is set to zero. Under these circumstances, a ther- 
mal gradient arises across the junctions and a station- 
ary heat current Jtot will flow from Si to S2, which are 
in steady-state thermal equilibrium. '^^H^^ As it was ar- 
gued in Ref. HH Jtot results from the sum of two terms, 
Jtot = Jqp{Ts^,Ts^) - Jint{Ts^,Ts^,ipa,^h,^c)- Here, 
Jqp is the heat current carried by quasiparticles that de- 
pends only on the temperatures of both superconduc- 
tors. On the other hand, Jint depends on the temper- 
ature of both superconductors and on the macroscopic 
phase difference across each junctioiP^ as well. This 
phase-dependent term is peculiar of weakly-coupled su- 
perconductors, and arises as a consequence of the inter- 
play between Cooper pairs and quasiparticles on tunnel- 
ing events through Josephson junctions. The Cooper pair 
condensate carries no entropy therefore leading to a zero 
contribution to the heat current .1^^ 
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In the following we shall concentrate on the phase- 
dependent heat current only. Taking into account that 
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FIG. 2. Top panels: phase-dependent component of the heat 
current Jint as a function of the magnetic flux in the main 
ring $1 plotted for different values of the control flux $2 for a 
double- loop heat interferometer with n = 1 and r2 = 1, i.e., 
Ra = Rb = Rc- The inset shows the case for r2 = 0. Bottom 
panels: transfer function 7j = dJint/d^i vs. $1 plotted for 
the same values of $2 as in the top panels. As it can be seen 
in panels (a) and (c), the amplitude of the oscillation 6 Jint 
remains constant for < $2 < ^o/3 decreasing for higher 
values of $2 [see panels (b) and (d)]. 



we are dealing with a three-junction circuit, Jint reads 



E 

=a,bA 



Jint(Tsi,Ts2) cos ifi 



(1) 

We emphasize that, whereas charge current de- 
pends on the sine of the phase difference across 
a Josephson junction,^ heat current depends 
on the cosine of cpi. In Eq. ([l]), J*^^ = 
^ J,^ desMrie, TsJM^ie, Ts, ) [ f{Ts,) - f{Ts, )] 
where Mj{e,Ts^) = A,{Ts^)/ ^e^ - A,{Ts^)^e[e^ - 
Aj{Ts^)% f{Ts^) = tanh(s/2feBT5,), Aj{Ts^) is the 
temperature-dependent energy gap of superconductor 
Sj^ Q{x) is the Heaviside step function, is the Boltz- 
mann constant and e is the electron charge. For the 
sake of completeness we also provide with the expression 
for the heat current carried by quasiparticles,^ 

Jqp{TSr^Ts^) = T.i=aAcJqp^^S^^Ts^)^ ^^^^^ 

Here ^j{s,Ts^) = \s\/ ^s^ - A,{Ts^ye[s' - A,{Ts^)^] 
is the BCS normalized density of states in Sj. 

The explicit dependence of Jint on the applied mag- 
netic fluxes and the circuit symmetry parameters, i.e., 
on the ratio between the normal-state resistances of the 
junctions, can be calculated analytically as follows. On 
the one hand we neglect the geometric inductance of each 
ring so that the fluxoid quantization on both of them im- 



(Pc + 27r^2/^o - 



2n7r 
2m7r, 
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Wb is the flux quantum, and 
The inclusion of the geometric 



where <l>o = 2.067 x 10 
n and m are integers, 
inductances of the loops does not contribute to the un- 
derstanding of the essential physics in our device, and 
makes the deduction of analytical expressions impossi- 
ble. The conservation of the circulating supercurrent in 
both loops, on the other hand, imposes 

ijsm(pa = i^jsin(piy — ijsm(pc^ (3) 

where ij is the Josephson current of junction i 
that, according to the generalized Ambegaokar-Baratoff 
model, is proportional to R~^' In writing Eq. ^ we 
have established a given current sign convention (see yel- 
low arrows in Fig. |T]). We define now the circuit symme- 
try parameters as ri = ij/i^j = Jint/ Jint = ^b/Ra ^ ^1 



and r2 



Jint/ Jint = Rb/Rc- ^2 can vary be- 



tween < r2 < 1 since setting r2 > 1 is equivalent to 
exchange the roles of and Rc. Combining Eqs. ([2| 
and ([3| and using simple trigonometric relations one gets 
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where a = cos(27r^i/^o), P = cos(27r^2/^o) and 7 = 
cos[27r(<l>i - ^2)/^o]' The choice of the signs of Eqs. (4) 
obeys the requirement of minimizing the free energy (Ej) 
of the whole system. The latter can be written as Ej = 
3Ej^o-Ej^o T.i=a,b,c^J^^^^i^ where Ej^o = ^o/27r.^^ By 
inserting Eqs. (4) into the previous expression it can be 
easily seen that the minimum of Ej is obtained for the 
solutions with positive signs. These, inserted into Eq. 
([1]), yield finally the following expression for Jint^ 



Jint = jL\/l + r? + r2 + 2ria + 2r2/3 + 2rir27. (5) 

We note that, if we set r 2 = 0, i.e., Rc ^ cxo, one recovers 
the expression corresponding to the single-loop heat in- 
terferometer conceived in Ref. [TH Consequently, for the 
case in which ri = 1 and r2 = 0, i.e., Ra = Rb and Rc 
00, Eq. (|5| becomes e_qual to Jint = "^Jintl cos(7r<l>i/<l>o)| 



[see the inset in Fig. [2|b)]. This function is maximized 
for <l>i = k^o and minimized for $1 = /c$o/2, k being an 
integer. The amplitude of the oscillation 5 Jint ^ defined 
as the difference between the maximum and minimum 
value of Jint, is given by 2Jf^^ in this case. 

We now turn our attention towards the double-loop 
heat interferometer discussed here. The straight est 
choice is ri = r2 = 1, i.e., Ra = Rb = Rc- Although 
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FIG. 3. Density plots showing J int as a function of the mag- 
netic flux on the main and the control loop, $i and $2, respec- 
tively. Three representative cases (denoted #1, #2 and #3) 
have been considered and are schematized on the top part of 
each plot. The maximum and minimum values of J int change 
perceptibly when reducing the symmetry of the double-loop. 



being the most simple configuration, it enables us to in- 
fer most of the characteristics of our thermal interfer- 
ometer. We can distinguish between two regimes, one 
defined by < <l>2 < *^*o/3 [Fig. |2|a)] and a second one 
covered by <l>o/3 < ^2 < ^o/2 pig. [2|b)]. If no mag- 
netic flux is applied to the control loop, i.e., $2 = 0, 
the phase-dependent term of heat current is given by 
Jint = -^fnt\/l + 8cos2(7r<l>i/<l>o). As showu in Fig. jija), 
this function exhibits the same expected <l>o-perioalcity 
and amplitude of the oscillation as that of the limit case 
in which ri = 1 and r2 = [see the inset of Fig. 2|b)] 
but the appearance of the oscillation is drastically differ- 
ent. In addition, notice that the minimum of J^nt does 
not go to zero anymore, i.e., there is no possibility of 
suppressing the phase-dependent component of heat cur- 
rent. As we increase the amplitude of the control flux ^2 
the mean value and the shape of the curves continue to 
evolve whereas 6 Jint holds unchanged. Furthermore, the 
curves turn out to be shifted horizontally. 

At <l>2 = <^>o/3 a noticeable phenomenon takes place. 
Under these circumstances Eq. ([5| takes the form Jint = 
2Jf^J cos[7r(<l>i/<l>o - 1/6)]|. This is to say, apart from 
a small shift equal to ^0/6, one recovers the same de- 
pendence on <l>i obtained for the symmetric single-loop 
heat interferometer. If we continue increasing the con- 
trol flux we enter in a new regime. Under these cir- 
cumstances it can be shown that SJint decreases linearly 
with ^2 whereas the mean value of Jint remains con- 
stant. Furthermore, at ^2 = ^o/2, the modulation dis- 
appears completely and Jint becomes independent of $1, 
i.e., Jint = Jint' The aforementioned characteristics are 
emphasized in Figs. |2|c) and |2|d) where we plot the 
transfer function, 7j = dJint/d^i^ for both regimes. 

It is worthwhile now to analyze what happens by re- 




FIG. 4. (a) Phase diagram for the double- loop heat inter- 
ferometer on the ri — r2 space. Within the white region the 
amplitude of the oscillation SJint increases linearly with ri. 
On the diagonally stripped region the maximum value of 6 Jint 
is obtained independently of ri and r2. Finally, on the grid 
region, SJint decreases with r2. Additionally, the two dashed 
lines define a triangle within which the interferometer can be 
fully balanced, (b) Jint vs. $1 for three selected points (#1, 
#3 and #4) contained within the aforementioned triangle and 
sketched on the right. In the three cases, the dashed line cor- 
responds to $2 = whereas the solid line corresponds to the 
value of $2 that provides the fully suppression of Jint- 



ducing the symmetry of our device. In Fig. [3] we show 
the density plots of Jint vs. ^1 and ^2 for three repre- 
sentative cases, including the symmetric double-loop that 
we have analyzed previously. In general, the maximum 
of Jint is always reduced for the cases in which one re- 
sistance is different from the others. Let us analyze in 
more detail what happens with SJint- If ^1 < 1 , i-e., 
Ra ^ Rb^ we find that the maximum amplitude of oscil- 
lation at fixed ^2 is given by SJint = "^Jint"^^ f^^ whatever 
value of r2, i.e., for whatever value of Rc- This is to say, 
SJint is independent of the degree of asymmetry of the 
control loop and increases linearly by increasing the sym- 
metry on the main loop. On the other hand, if ri > 1 
the dependence of SJint at fixed ^2 is more complicated. 
For ri — r2 < 1, 5 Jint is independent from the degree 
of symmetry in both loops reaching it maximum value, 
that, in this case, is given by SJint = "^Jint- Finally, for 
Ti — r2 > 1, SJint decreases with r2- These regimes are 
summarized in Fig. [4]^ a). 

Let us ascertain now whether it is always possible to 
suppress completely Jint- The minimum of Jint can 
be easily determined from Eq. ([5|. By requesting 
dJint/d^i = and d Jint/ 9^2 — 0, and by summing 
the resulting equations we get that ri sin(27r<l>i/<l>o) + 
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r2 sin(27r^2/^o) = 0- This equation gives us the family 
of values of $i and $2 that maximize or minimize Jint 
for any given ri and r2. For instance, if $1 = <^>o/2, 
there are only two solutions, <l>2 = and ^2 = *^o/2, 
that satisfy this condition and, in addition, correspond 
to a minimum. Inserting these solutions into Eq. ([5| and 
imposing Jint = we obtain ri — r2 = 1 and ri + r2 = 1, 
respectively. These equations define two straight lines in 
the ri — r2 space plotted in Fig. |4] Within this area 
there exists at least one value of ^2 that enables us to 
write Jint as a function of | cos [7r(^ 1/^0 ~ ^)]| where is 
a shift in <l>i. The aforementioned conditions take a more 
eloquent form when expressed in terms of the Josephson 
critical currents of each junction, giving 

ia - ic < ib < ia ^ ic (6) 

Unlikely to a single- loop heat interferometer,^^ even a 
quite asymmetric double-loop structure for which in- 
equality ([6| holds, offers the possibility of suppressing 
completely Jint through an appropriate choice of ^2- 

The three cases plotted in Fig. [3] satisfy the aforemen- 
tioned conditions and shall therefore be useful to illus- 
trate this behavior. In Fig. ^h) we plot Jint{^i) for 
^2=0 (dashed lines) and for the corresponding value of 
^2 that provides the fully suppression of Jint (solid lines). 
The curves corresponding to the symmetric double-loop 
have already been plotted in Fig. [2|a) and[2|b), and are 
therefore not shown. Notably, when setting ri = 0.5 and 
r2 = 0.5, i.e., Ra = Rc = '^Rb^ Jint cancels precisely 
when ^2 = *^o/2 at = <^>o/2 ± k^o since this case 
satisfies exactly the condition ri + r2 = 1. Notice that 
the total amplitude of the oscillation is reduced by one 
half with respect to the previous example. Let us finally 
consider a last illustrative case belonging to the region 
defined by 1 < ri < r2 + 1, i.e, the diagonally stripped 
region in Fig. Qa). If ri = 1.25 and r2 = 0.75, i.e., 
Ra = 0.6Rc and 2i?5 = 0.75i?c, although corresponding 
to a substantially asymmetric interferometer, it should 
be possible to suppress completely Jint while conserv- 
ing the maximum amplitude of oscillation. As we can 
see in the bottom panel of Fig. ^h) this is exactly the 
case. Setting ^2 = *^*o/4, the phase-dependent compo- 
nent of heat transport cancels at $1 = 3<l>o/5±/c^o with 

^Jint = '^Jinf 

We shall finally dedicate a few words to some potential 
applications and practical aspects related to the fabrica- 
tion of the heat interferometer proposed here. Our struc- 
ture can be integrated within, not only superconducting 
elements, but also hybrid mesosocopic circuits composed 
of, e.g., normal metals, two dimensional electron gases 
and semiconductor nanowires as well. A precise control 
of the amount of heat flowing through such circuits is of 
crucial importance.'!'^ Temperature determines, for in- 
stance, the phase transition in superconductors,!^ the 
amount of heat exchanged between electron and lat- 
tice phonons,'!' the energy level occupation in quantum 
systems,^ or the critical current flowing through Joseph- 
son junctions.^ Mastering the heat current in supercon- 



ducting circuits w ould enable the in-situ fine tuning of 
radiation detectors.^^ Controlling the temperature of a 
two-level quantum system can eventually have influence 
on its decoherence time or contribute to its initialization 
in quantum computing architectures.^ Furthermore, fully 
tunable Josephson junctions of different kinds can be en- 
visioned. In such devices, the direct relation between 
the electronic temperature and the critical current can 
be exploited to modulate the latter via the application of 
a magnetic flux. Unlike the usual voltage- controlled hot- 
electron Josephson transistors,'^^H^the principle of oper- 
ation proposed here can lead to magnetic flux- controlled 
thermal Josephson transistors. 

A single- loop double-junction superconducting ther- 
mal device connected to two normal-metal electrodes 
has indeed recently been realized. The existence of 
the phase-dependent component of heat current Jint has 
therefore been proved experimentally. Our double-loop 
heat interferometer scheme provides further advantages 
whereas it does not imply extra difficulties from the 
point of view of the fabrication. Such structure can be 
easily fabricated by standard electron-beam lithography 
and shadow mask evaporation of superconducting met- 
als, e.g., aluminum. Aluminum oxide for the Joseph- 
son barriers can be used. The realization of analogous 
Josephson devices composed of two or more loops with 
independent magnetic flux controls operat ing a s charge 
interferometers has been recently reported] ^ 1 ^^ ^ proving 
the feasibility of this structure. 

To conclude, we have provided with all the informa- 
tions required for designing a fully-balanced heat inter- 
ferometer. Such a device allows, on the one hand, to 
modify the form and phase shift of the phase-dependent 
heat current by tuning the control flux. This would en- 
able the user to choose a convenient point of operation 
within the available flux-to-heat current transfer charac- 
teristic. On the other hand, we have demonstrated that 
if condition (|6| holds, a fully-balanced interferometer is 
obtained, meaning that the phase-dependent part of the 
heat current can be completely annihilated. 
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